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Abstract 
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$— i In this note we prove the exact formula for the expected number of moves to solve several variants 

Mh of the Tower of Hanoi puzzle with 3 pegs and n disks, when each move is chosen uniformly from 

^ the set of all valid moves. 

CO 

1 Introduction 

We define two variants of the Tower of Hanoi puzzle with the following starting and final positions: 



O 

u 

Puzzle 1. The starting position is uniformly chosen from the set of all 3 n positions. The final 
position is with all disks on the same peg (no matter which one). 



Puzzle 2. The starting position is with all disks on the first peg. The final position is with all disks 
on the third peg. 



Let E\(n) and E2(n) be the expected number of random moves required to solve Puzzle 1 and 
qq 2 respectively. 

Puzzle 1 was posed by David G. Poole who submitted a sequence Ei(n) for n up to 5 to the 
OEIS [ ]. Later Henry Bottomley conjectured the following formula for Ei(n)\ 

m E n n ) = - A • W 



Puzzle 2 was posed by the second author [ ] who also submitted to the OEIS [! ] values of ^(n) 
for n up to 4 but did not conjecture a general formula for E2(n). 



In this note we will prove the formula (1) as well as the following formula for E2(n): 



(3 n -l)(5 n -3 n ) 
2-3 r 



E,(n) = '" -■'£-,•" P> 



originally found by the first author. 

Let Dk be the disk of size k such that in n-disk puzzle D\ and D n refer to the smallest and 
largest disks respectively. Similarly, let D™, where k > m, be the set of all disks of sizes from m 
to k inclusively. 
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2 Auxiliary Puzzles 

We start by defining two auxiliary puzzles: 

Puzzle 3. The starting position is with all disks on the first peg. The final position is with all disks 
on the same peg. At least one move is required. 

Let Es(n) be the expected number of moves required to solve this puzzle, while pi(n) and 
p 2 (n) denote the probability that a final configuration has all disks on the first and second peg, 
respectively. 

From the symmetry, it is clear that the probability that a final configuration has all disks on 
the third peg is also p 2 (n), so pi(n) + 2p 2 (n) — 1. 

Puzzle 4. The starting position is with the largest disk on the second peg and the other disks on 
the first peg. The final position is with all disks on the same peg. 

Let £4(71) be the expected number of moves required to solve this puzzle, while gi(n), q 2 (n), 
and qs(n) denote the probability that a final configuration has all disks on the first, second, and 
third peg, respectively. 

The relationship between Puzzles 3 and 4 and Puzzles 1 and 2 is given by the following theorems. 

Theorem 1. 

£i(n) = £i(n-l) + -£ 4 (n) 

Proof. Note that D n cannot move unless Dj l _ 1 are on the same peg. In Puzzle 1, the expected 
number of steps required to arrive at such position is E\(n — 1). Moreover, since the starting 
position is uniformly chosen, Dj l _ 1 will be on any peg with the equal probability 1/3. In particular, 
if it is the same peg where D n resides, then the puzzle is solved. Otherwise, we can view the 
remaining moves as solving Puzzle 4, implying that 

Ei{n) = -£i(n - 1) + -(£i(n - 1) + £ 4 (n)) = £i(n - 1) + -£ 4 (n). 



□ 



Theorem 2. 

E 2 {n) 



E 3 (n) 
p 2 {n) 



Proof. Note that D n cannot move unless Dj l _ 1 are on the same peg. 

In solution of Puzzle 2 with random moves, first time all disks will appear on the same peg 
after £3(71) moves on average. This peg will be the first one with the probability pi(n), the 
second one with the probability p 2 (n), or the third one with the probability P2{p). In the last case 
the puzzle is solved, while in the first two cases we basically obtain a new Puzzle 2. Therefore, 
E 2 (n) = E 3 (n) + (pi(ra) +p 2 (n))E 2 (n), implying that E 2 {n) = gg. □ 

Therefore, to find explicit formulae for Ei(n) and E 2 (n), it is enough to find those for E 3 (n) 
and £4(71) along with the probability p 2 (n). 



3 Derivation 

Lemma 3. The following equalities hold: 
(i) E 3 (n) = E 3 (n - 1) + 2p 2 {n - l)E 4 (n) 

(ii) pi(n) = pi(n - 1) + 2p 2 {n - l)q 2 (ri) 

(Hi) p 2 {n) =p 2 (n- l)gi(n) + p 2 (n- l)© (™) 

^ £7 4 (n) = \ + £ 3 (^ - 1) + (pi(n - 1) + p 2 (n - l))E 4 (n) 

(v) qi{n) =p 1 (n- l)gi(n) + p 2 (n- 1)?3 H 

(^ g 2 (n) = f(p2(rc- 1) + (Pi(^~ 1) +P2(n- l))q 2 (n)) + |(pi(n- l)qs(n) + p 2 (n - l)gi(ra)) 

fra,) g 3 (n) = |(pi(n- 1)93 H + p 2 (n- 1)<7iH) + \({pi(n - 1) + p 2 (n- l))g 2 (n) + p 2 (n - 1)) 

Proof. Consider Puzzle 3. Note that D n cannot move unless D\_ 1 are on the same peg. Therefore, 
we can focus only on D\_ 1 until they all come to the same peg, that will happen (on average) 
after E 3 (n — 1) moves. When D\_ 1 are on the same peg, it is the first one (where D n is) with the 
probability pi(n — 1), in which case we have the final position with all disks on the first peg. If this 
is not the case (with the probability 1 — pi(n — 1) = 2p 2 {n — 1)), we have D n on the first peg and 
D\_ 1 on a different peg (equally likely on the second or the third one). Therefore, the remaining 
moves can be considered as Puzzle 4 and, hence, the expected number of remaining moves in this 
case is E±(n). Therefore, the expected total number of moves to solve Puzzle 3 is 

Pl (n - 1) • E 3 (n - 1) + 2p 2 (n - 1) • (E 3 (n - 1) + E 4 (n)) = E 3 (n - 1) + 2p 2 {n - l)E 4 (n) 

which proves formula (i). 

From the above it is also easy to see that in the final position all disks will be at the first peg 
with the probability pi(n — 1) + 2p 2 (n — l)q 2 (n) and at the second peg or third peg with the same 
probability p 2 (n — l)q\{n) +p 2 (n — l)q 3 (n), that proves the formulae (ii) and (hi). 

Now, consider Puzzle 4. Moves in this puzzle can be split into two or three stages as follows. 
In Stage 1 only D n is moving (between the second and third pegs), Stage 2 starts with a move of 
D\ (from the top of the first peg) and ends when D\_ 1 are on the same peg. If this is not the final 
position, the remaining moves are viewed as Stage 3. Let us analyze these stages. 

It is clear that the expected number of moves in Stage 1 is 1 /2 while at the end of this stage 
we have D n on the second peg with the probability 3 / 4 an d on the third peg with the probability 
1 /a. The expected number of steps in Stage 2 is simply E 3 (n — 1) and at the end D\_ 1 are on the 
first peg with the probability p\{n — 1) and on the second or third pegs with the equal probability 
p 2 {n — 1). Therefore, with the probability p 2 (n — 1) we are at the final position (no matter where 
D n is left after Stage 1) and with the probability 1 — p 2 (n — 1) — pi(n— 1) J rp 2 (n — 1) we continue 
Stage 3. Stage 3 can be simply viewed as a new Puzzle 4 with the expected number of moves E 4 (n). 
The above analysis proves formulae (iv)-(vii). □ 

Note that the formula for E 3 (n) easily follows directly from (i) and (iv). Namely, the formula 
(iv) can be rewritten as p 2 (n — l)E 4 (n) — | + E 3 (n — 1) and substituted into (i), resulting in the 
following recurrent formula 

E 3 (n) = E 3 (n - 1) + 2p 2 (n - l)E 4 (n) = 3E 3 (n - 1) + 1. 



Together with £3(1) = 1 it implies 



3 n - 1 3 n - ] 



2p 2 (n-l)' 

Let us focus on the recurrent equations (ii), (hi), (v), (vi), (vii) and solve them with respect 
to pi(n), ^2(^)5 9i(^) 5 #2(^)5 an d Q3(n). Solving Puzzles 1 and 2 for n — 2, we have the following 
initial conditions: 

Pl(2) = 5/ 8 , p2(2)=3/i 6 , 9i(2) = l/8, 92(2) =5/8, 93(2) = 1/4- 

Also solving Puzzle 1 for n = 1, we get pi(l) = and ^2(1) = V 2 - 

From (v) we have P2(n — l)qs (n) = #i(n)— pi(n — l)qi(n) = {l—p\{n—l))qi(n) — 2p2(n — l)qi(n) 
and since p2(n—l) > for all n > 2, we have #3(71) = 2gi(n), ^(^) = 1 — 9i(^) —Q3(n) = 1 — 3gi(n), 
and (8 — 3pi(n — l))#i(n) = 1 for all n>2. 

Using these relations, we simplify equation (ii) to 

pi(n) = pi(n - 1) + 2p 2 (n - 1)(1 - 3<?i(n)) = 1 - 6p 2 (^ - l)?i(^) 
= 1-3(1 -pi{n- l))qi(n) = 1 -3gi(n) + 3p 1 (n- l)gi(n) 
= 1 - 3gi(n) + 8gi(n) - 1 = 5gi(n). 

Combining the above equations, we have (8 — 15gi(n — l))#i(n) = 1, implying that qi(n) = 
f( n \0) (i.e., the n-th iteration of /(•)) where 

m ' 



8 -15a; 
It is easy to prove by induction that 



Therefore, 



?i 


5- 1 - 
\ n ) = — ^ 


. gn-. 

. 3 n 


L 


?2(n) 


= 1 -3qi(n) : 


2- 


5 n-l 


5 n 


_ 3 n' 


93 H = 


= 2gi(n) = 2- 


n-1 _ 

5 n - 


_ 3 n-l 

_ 3 n 


pi(n) = 


5 71 
= 5gi(n) = 


'-5 
5 n - 


3 n 


7)0^79,^ 


__ 1 -pi(ra) 


3 


n-1 



5 n _ 3 n ' 



and finally 



^4(n) = ^(5- 1 -3"- 1 ) 



Now we are ready to prove formulae (1) and (2). Theorem 1 together with Ei(0) — implies 

5 n - 2 • 3 n + 1 
3 ~ 6 = 4 • 

k=l k=l 

Theorem 2 implies 



E^n) = Y, E i( k ) ~Ei(k- 1) = Y,^' 1 ~ 3fc_1 

fe=i 

£ 3 (n) (3 n -l)(5 n -3 r 



£ 2 (n) 



P2(n) 2 • 3 



re-l 



4 One More Puzzle 

Let us consider yet another puzzle somewhat similar to both Puzzles 1 and 2: 

Puzzle 5. The starting position is uniformly chosen from the set of all 3 n positions. The final 
position is with all disks on the first peg. 

Solving this Puzzle can be viewed as first solving Puzzle 1 and if it does not result in all disks 
on the first peg (that happens with the probability 2 /3), continue solving it as Puzzle 2. Therefore, 
the expected number of moves in Puzzle 5 is: 



n 



„ , . „ , . 2^ /x 5 n -2-3 n + l (3 n -l)(5 n -3 n ) 5 n+1 - 2 • 3 n+1 + 5 /5 
E b (n) = E^n) + -E 2 (n) = + — = 

5 Analysis of Puzzle 2 via Networks of Electrical Resistors 

We now present an alternative analysis of Puzzle 2 that couples a theorem about expected commute 
times of random walks on graphs with the delta-to-wye transformation originally introduced to aid 
in the analysis of three-phase AC systems for electrical power distribution. Building on a monograph 
by P. G. Doyle and J. L. Snell [ ], A. K. Chandra et al. [ >] proved the following theorem: 

Theorem 4 (The Mean Commute Theorem). The expected number of steps in a cyclic random 
walk on an undirected graph that starts from any vertex A, visits vertex B, and then returns to A 
equals 2i7iRab, where m is the number of edges in the graph and Rab is the electrical resistance 
between vertices A and B when a 1-ohm resistor is inserted in every edge of the graph. 

A Hanoi Tower with n = 1 may be represented by a graph with three nodes, which we shall 
denote by A, £>, and C, and three egdes AB,BC, and AC. A, B and C correspond, respectively, 
to the three possible states, U D\ is on peg 1", "D\ is on peg 2", and "D\ is on peg 3". Edge 
AC's presence represents the fact that it is possible to reach C from A, or A from C, in one move; 
the other two edges have analogous interpretations. Now insert a 1 ohm resistor in each of these 
edges. There are two parallel paths between A and C, the direct path along edge AC which has a 
resistance of 1 ohm, and the indirect path from A to B to C that has a resistance of 2 ohms, so the 
overall resistance from A to C is 2 /3 ohm. Since there are 3 edges in the graph, the mean commute 
time from A to C and back is 2 • 3 • 2 /s = 4. By symmetry, on average half of the time is spent 
going from A to C and the other half returning from C to A. Accordingly, the mean time it takes 
a randomly moving Hanoi Tower with n = 1 to reach peg 3 starting from peg 1 is 4 /2 = 2, which 
agrees with our formula for E^(^) when n = 1. 

We proceed to iterate this approach in order to derive the formula for E2(ri) for general n. The 
key is to use the classical delta-to-wye transformation of electrical network theory [ ]. A "delta" 
is a triangle with with vertices A, B and C that has resistances tab hi edge AB, tac m edge AC, 
and tbc m edge BC, as depicted in Fig. la. The corresponding "wye," shown in Fig. lb, has the 
same three nodes - A,B, and C - plus a fourth node, x, and three edges Ax, Bx and Cx that contain 
respective resistances r a ^r\) and r c . It is straightforward to verify that, if 

r — ^AB^AC r — ^AB^BC an( J r — ^AC^BC (0\ 

a ^AB+^AC+^BC ' b ^AB+^AC+^BC C ^AB+TAC+^BC ' ^ ^ 

then the net resistance Rab between nodes A and B will be the same in Fig. lb as it is in Fig. la, 
and likewise for the net resistances Rac between nodes A and C and Rbc between nodes B and 
C. 
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Figure 1: 

We shall need to consider only the special case tab = ^AC — ^BC — R, m which r a — r^ — r c — 
R /3. In particular, when R = 1, we have r\ — r^ = V 3 ? so Fig. lb yields Rac — 2 / 3 ? the same result 
we obtained before by considering the two parallel paths from A to C in Fig. la. 

The graph for the state transitions of the Hanoi Tower for n — 2 is the Sierpinski gasket, or 
fractal, shown in Fig. 2. It is obtained by thrice replicating the graph for n — 1 and then connecting 
each of the resulting three pairs of non-extreme vertices by a bridging link. This leaves only three 
extreme vertices which we continue to label by A, B and C and which correspond to both disks 
being, respectively, on peg 1, on peg 2, and on peg 3. The other nodes, labeled a,/3, 7,5, A and 
/x, correspond respectively to the following ordered triples of contents of peg 1, peg 2 and peg 3: 
a = (Z> 2 , 0, Di), /3 = (Z>2, D u 0), 7 = (0, D 2 , Z?i), 5 = (0, D U D 2 ), A = (D u D 2 , 0), M = (A, 0, ^2), 

Theorem 5. Delta-to-Wye Induction 

The Sierpinski gasket state diagram for an n-disk Hanoi Tower with a unit resistance in each of 
its branches can be converted, for purposes of determining the resistance between any two of the 
extreme vertices A, B and C, into a simple "wye" in which A, B and C each are connected to a 
centerpoint x, by links that each contain a common amount of resistance we shall denote by Rn/2. 

Proof. We now prove Theorem 5 by induction. We have already shown that it is true for n — 1, 
the value of ^1/2 being !/3 ohm. We assume the theorem is true for n equal to the positive integer 
k and proceed to show that it them must also true for n = k + 1. 

The Sierpinski gasket graph for a Hanoi Tower with k + 1 disks is produced by generating three 
replicas of that for a Hanoi Tower with k disks that possess respective extreme vertices {A\, £>i, C\\. 
{A 2 , £>2, C 2 } and {A3, S3, C3} and then adding three bridging links, as shown in Fig. 3, one between 
B\ and A 2 , another between C\ and A3, and the third between C 2 and S3. The extreme vertices in 
the resulting graph are A\^B 2 and C3, i.e. the only A^'s, B^s and Q's to which none of the three 



A 







Figure 2: 

bridging links is incident. When one seeks to apply Theorem 4, a unit resistance is inserted in each 

of the three bridging links, just as is the case for every other link in the graph. By the induction 

hypothesis, for purposes of computing the resulting resistances between any two of the vertices 

Ai, £>i, Ci, A2, £>2, C2, A3, S3, and C3, the remainder of the resistive graph for n = k + 1 consists of 

three wye's, each comprised of three links each of which have resistance R k/2, as shown in Fig. 4. 

Denoting the centers of these three wye's by x, y, and z, we see that triange xyz is a delta each 

edge of which has resistance R k/2 + 1 + R k/2 = R^ + l. Applying the delta-to-wye transformation to 

this delta network results in the wye network of Fig. 5. Each of the links in this wye has the same 

resistance, namely 

R h±1 = R, R^l = 5_ 1 

2 2 3 6 3 y J 

Theorem 5 is proved. □ 

In the process of proving Theorem 5, we have obtained the first order linear diffrence equation 

Rk+i = -~Rk + o- ( 5 ) 

From this equation and the boundary condition R\ — 2 /3, we obtain the key result, 

R n=(l) "I- (6) 

The number, m n , of edges in the Sierpinski gasket graph for an n-disk Hanoi Tower is 

^ = (3--3M + 3.2 = 3 (3 „_ 1) 



A=A< 
A 



/ 9&i 



Az 



/ 



f 



\ 




0>»*pnsi*\ W SG^s o** Him 



Figure 3: 

From the Mean Commute Time theorem and the symmetry of a random walks from A to C and 
from C to A, it follows that the mean number of steps it takes a randomly moving n-disk Hanoi 
Tower to transfer all its disks from peg 1 to peg 3 equals 

(3 n -l)(5 n -3^ 

This concludes our alternative derivation of formula (2) for the mean number E2(n) of moves in 
Puzzle 2. 



m n R A c = m n R n = 1 . (7) 



6 Historical Commentary 

The French mathematician, Edouard Lucas, created the Tower of Hanoi puzzle in 1883 [ ]. Appar- 
ently, he simultaneously created the following legend [ ]: 

"Buddhist monks somewhere in Asia are moving 64 heavy gold rings from peg 1 to peg 
3. When they finish, the world will come to an end!" 



'v 2 ? Zl. \ f Rk/2 
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Fig. # A network of resistors apropos of a 
Hanoi Tower with Jc+t disks 



Figure 4: 

That Hanoi is located in what was then French Indo- China perhaps explains why a Frenchman saw 
fit to include Hanoi in the name of his puzzle. However, the legend never placed the monks and 
their tower explicitly in Hanoi or its immediate environs. Lucass Hanoi Tower puzzle became an 
international sensation (think Loyds 15 puzzle and Rubiks cube); the legend was used to bolster 
sales. Still a popular and beloved toy, the Hanoi Tower now also can be accessed over the Internet 
as a computer applet. 

Since it is often intimated that monks possess superhuman abilities, maybe they can move 
large gold objects rapidly. Perhaps they can make a move a microsecond, maybe even a move a 
nanosecond. Considering that for n — 64 the minimum number of moves required to solve Puzzle 
2 is "only" 2 64 1 = 18, 446, 744, 073, 709, 551, 615, planet Earth may expire any day now. This is in 
part what motivated suggesting adoption of a randomly moving Hanoi Tower [ ]. Formula (7) shows 
that replacing the minimum-moves strategy with a random walk forestalls the end of the world by 
a factor of roughly (|) > 2.9 x 10 25 on average. Although this is reassuring, it nonetheless would 
be further comforting to know that the coefficient of variation of the random number of steps in 
Puzzle 2 with n — 64 is small, i.e., that its standard deviation is many times smaller than its mean. 
Exact determination of said coefficient of variation is an open problem that we may address in 
future research. An extensive bibliography of some 370 mathematical articles concerning the Hanoi 
Tower puzzle and variations thereon has been complied by Paul Stockmeyer [ ]. 



yz 










H^uerV 



A 




-$Vi~&L 



5 R j. J- - P 






Figure 5: 

References 

[1] T. Berger. Lucas, Sierpinski, Markov, Shannon and The End of the World. In Invited presen- 
tation in the David Slepian Memorial Session, Information Theory and Its Applications (ITA 
2008), UCSD, La Jolla, CA, January 2008. 

[2] Toby Berger. Sequence A134939 in N. J. A. Sloane (Ed.), The On-Line Encyclopedia of Integer 
Sequences. 2008. http://oeis.org/A134939. 

[3] A. Chandra, P. Raghavan, W. L. Ruzzo, R. Smolensky, and P. Tiwari. The electrical resistance 
of a graph captures its commute and cover times. Seattle, May 1989. ACM Press. 

[4] P. G. Doyle and J. L. Snell. Random Walks and Electrical Networks. The Mathematical 
Association of America, Washington, D.C., 1984. 

[5] A. E. Kennelly. Equivalence of triangles and stars in conducting networks. Electrical World 
and Engineer, 34:413-414, 1899. 

[6] E. Lucas. Recreations Mathematiques, volume III. Gauthier-Villiars, 1893. Reprinted by Albert 
Blanchard, Paris. 

[7] David G. Poole. Sequence A007970 in N. J. A. Sloane (Ed.), The On-Line Encyclopedia of 
Integer Sequences. 2008. http://oeis.org/A007798. 



[8] 



R. L. Ripley. The New Believe It or Not Book 2nd Series. Simon and Schuster, 1931. 



[9] P. K. Stockmeyer. The Tower of Hanoi: A Bibliography. Version 2.2. 2005. http://www.es. 
wm . edu/~pkstoc/biblio2 . pdf . 



